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ABSTRACT 

We present two transits of the hot- Jupiter exoplanet XO-2b using the Gran Tele- 
scopio Canarias (GTC). The time series observations were performed using long-slit 
spectroscopy of XO-2 and a nearby reference star with the OSIRIS instrument, en- 
abling differential specrophotometric transit lightcurves capable of measuring the exo- 
planet 's transmission spectrum. Two optical low- resolution grisms were used to cover 
the optical wavelength range from 3800 to 9300 A. We find that sub-mmag level slit 
losses between the target and reference star prevent full optical transmission spectra 
from being constructed, limiting our analysis to differential absorption depths over 
~ 1000 A regions. Wider long slits or multi-object grism spectroscopy with wide 
masks will likely prove effective in minimising the observed slit-loss trends. During 
both transits, we detect significant absorption in the planetary atmosphere of XO- 
2b using a 50 A bandpass centred on the Nai doublet, with absorption depths of 
A(R p i/R±) 2 = 0.049±0.017% using the R500R grism and 0.047±0.011% using the 
R500B grism (combined 5.2-ct significance from both transits). The sodium feature 
is unresolved in our low-resolution spectra, with detailed modelling also likely ruling 
out significant line- wing absorption over an ^800 A region surrounding the doublet. 
Combined with narrowband photometric measurements, XO-2b is the first hot Jupiter 
with evidence for both sodium and potassium present in the planet's atmosphere. 
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1 INTRODUCTION 

Using transiting exoplanets to obtain transmission spec- 
troscopy is now a proven method to probe planetary at- 
mospheres. To date, most of the observed atoms and 
molecules identified with this technique have been for the 



* Based on observations made with the Gran Telescopio Canarias 
(GTC), installed in the Spanish Observatorio del Roque de los 
Muchachos of the Instituto de Astrofsica de Canarias, in the is- 
land of La Palma, and part of the large ESO program 182.C-2018. 
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two most favourable cases, the hot Jupiters HD 209458b 
and HD 189733b. The detections include species in the up- 
per escaping atmosphere such as hydrogen, oxygen, carbon, 
and silicon \ Vidal-Madjar et al]|2003| |2004| |Lecavelier Des 



Etangs et al.|2010||Linsky et al.|2010|| Jensen et al.|2012[ ). In 

the lower atmosphere, sodium has been detected in both of 
these planets (|Charbonneau et al.|2 002 Redfield et al.|20"08 



Snellen et al. |2008 1 with the line profile also able to place 
constrains on the temperature-pressure profile (|Sing et al.| 



|2008||Vidal-Madjar et al.|2011b|a||Huitson et al.|2012[ ). Al- 

kali metals in other hot Jupiters have been detected as well, 



© 2012 RAS 



2 D. K. Sing 



with evidence for sodium in Wasp-17b (Wood et al. 20111 
and potassium in XO-2b ( Sing et al.|2011a l. 

Contrary to clear hot-Jupiter atmospheric models, 
HD 189733b shows ample evidence for a dominant haze 
which extends from the optical ( Pont et al.|2 008 Sing et al 



2011b[) into the near-IR ( |Sing et al.|2009||Desert et al.|2009 



2011; Gib son et al.|2012| ). The signature of the haze is char 



acteristic of Rayleigh scattering by small particles ( Lecave- 
lier des Etangs et al.|2008 1 , with silicate grains a possibility 



as predicted by some cloud mo dels ([Marley et al.||1999| |Su-| 
darsky et al.|2000| [Seager et al.||2000| |Helling et al.|2008[ ). 

In the context of hot Jupiters, the alkali metals sodium 
and potassium are extremely important. Alkali metals are 
thought to be responsible for the low albedos of hot Jupiters, 
with the pressure broadened potassium and sodium line- 
wings efficiently absorbing the incoming optical stellar ra- 
diation (ISudarsky et all 120001 ISeager et all 120001). Most 



hot Jupiters do appear to have low albedos (Cowan & Agol 



2011 ), though there are exceptions (Demory et al. 2011 San 



terne et al|2011 1. As the alkali metals are thought to play 
such a dominant role in the planetary energy budget of hot 
Jupiters, dictating where the stellar radiation is deposited, 
placing observational constraints on these metals is an im- 
portant first step. 

Sodium and potassium can both be detected from 
space-based or ground-based spectra or photometric mea- 



surements (e.g. 


Charbonneau et al. 2002 Redfield et al. 


2008 Sing et al. 


2011al. The challenge of obtaining broad- 



band transmission spectroscopy from the ground is the re- 
quirement for very high photometric accuracy, along with 
good spectral resolution on both the target and reference 
stars. Observations toward this end have been made on large 



aperture telescopes with integral field spectroscopy ( Arribas 



et al. 


2006) 


2010 


2011). 



2006) and multi-object spectrographs (Bean et al. 
20111. Here, we report the use of long-slit time-series 
spectroscopy to take simultaneous measurements of a tran- 
sit event along with a nearby reference star, which is used 
for differential spectrophotometry. Such multi-object spec- 
troscopic techniques have the potential to deliver transmis- 
sion, emission, and albedo spectra of transiting extrasolar 
planets, making such observations competitive with space- 
based instruments. The larger telescope apertures available 
from the ground also lead to better time sampling of the 
short hour-long transit events and the ability to observe 
fainter targets, an important step in being able to charac- 
terise the large number of fainter transiting planets discov- 
ered from photometric surveys. The observations reported 
here are part of a large ESO programme on the GTC, per- 
forming a spectrophotometric optical survey of transiting 
hot Jupiters using narrowband imaging and long-slit spec- 
troscopy. In this paper, we describe the GTC XO-2b ob- 
servations in §2, present the analysis of the long-slit transit 
light curves in §3, discuss the results and compare to model 
atmospheres in §4, and conclude in §5. 



2 OBSERVATIONS 

We observed the exoplanet host star XO-2A with the GTC 
10.4 metre telescope installed in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astrofisica 
de Canarias on the island of La Palma, using the Optical 



System for Imaging and low Resolution Integrated Spec- 
troscopy (OSIRIS) instrument during two separate nights 
in service mode. 



2.1 GTC OSIRIS instrument setup 

The GTC OSIRIS instrument ( |Cepa et al.|2000| |2003[ ) con- 
sists of two Marconi CCD detectors, each with 2048x4096 
pixels and a total FOV of 7.8' x8. 5', giving a plate scale of 
0.127"/pixel. The CCD has a pixel well depth of ~100,000 
electrons and a 16 bit A/D converter, which saturates at 
65,536 counts. For our observations, we choose the 2x2 bin- 
ning mode with a full frame readout at a speed of 500kHz, 
which produces a readout overhead time of 17 seconds. The 
readout speed of 500kHz has a gain of 1.46 e~/ADU and a 
readout noise of 8 e~ . While 500kHz is the noisiest readout 
speed (100 and 200kHz are available), it allows for the most 
counts per image to be obtained due to a higher gain and a 
shorter readout time which provides a better sampling rate. 

We utilised the spectroscopic capabilities of OSIRIS, 
using a long slit of 8.67' in the spatial direction with a 5" 
width. We used the widest available slit at the time of the 
observations to minimise slit losses. We observed two transit 
events of XO-2b, observing one with the R500B and one with 
the R500R grism. The R500B grism covered a spectral range 
from 3750 to 8586 A, has a dispersion of about 3.7 A per 
binned pixel, and has a peak efficiency of 17% near 5000A. 
The R500R grism covered a spectral range from 5000 to 9300 
A, has a dispersion of about 4.7 A per binned pixel, and has 
a peak efficiency of 15% near 6000A. 

For the XO-2b observations, we placed both the target 
and a chosen nearby reference star, XO-2A's stellar compan- 
ion XO-2B, within the slit, orientating the spatial direction 
such that both stars were centred within the 5" slit (see 
Fig. [I}. The companions XO-2A and XO-2B are separated 
by only 31" and cataloged as each having the same spectral 
type (K0V) and have very similar magnitudes and colours, 
making XO-2B an ideal reference star. In the future, multi- 
object spectroscopy modes for OSIRIS are planned, which 
will provide an added benefit over long-slit spectroscopy as 
more reference stars will be obtainable and wider slit masks 
can be fabricated. In addition, wider long slits are planned 
for OSIRIS as well. 



2.2 Observing log 

XO-2A and XO-2B were observed keeping them simulta- 
neously on the slit during the XO-2b transit events of 23 
December 2011 and 23 February 2012. 

23 Dec. 2011: The seeing was stable during the observa- 
tions, with values ranging from 0.8" to 0.9" under photomet- 
ric conditions. The R500R grism was used with 10 second 
exposure times, with 770 exposures obtained between UT 
01:01 and UT 06:56, with airmass values ranging between 
1.07 and 1.54. 

23 Feb. 2012: During the night, the seeing ranged from 
0.8" to 2.0" under photometric conditions. The R500B grism 
was used with 22 second exposure times, with 519 exposures 
obtained between UT 20:32 and UT 02:15, with airmass val- 
ues ranging between 1.07 and 1.4. 

For the two nights, exposure times were not modified 
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Figure 1. GTC Osiris CCD frame of XO-2A (right) and XO-2B 
(left) long-slit spectra from 23 Feb 2012 taken with the R500B 
grism and a 5" wide slit. The spectra cover wavelength ranges 
from 3750 A (bottom) to 8587 A (top), with background sky 
emission lines seen in the cross-dispersion direction. 



during the sequence, maintaining the peak counts between 
25000-40000 ADUs for both the reference and target stars, 
well within the linear regime of the detectors. When needed, 
defocusing was used to avoid saturation in improved seeing 
conditions, with the full width at half maximum (FWHM) of 
the point spread function (PSF) typically kept about ~2.5 
times smaller than the slit width. At the expense of higher 
spectral resolution, adopting a modest defocus also helps to 
minimise changes to the spectral resolution. For the 23 Feb. 
observations, 97% of the spectra have a FWHM within 25% 
of the average, despite a factor of more than two variation in 
the reported seeing conditions. The wavelength calibration 
was determined from HgAr, Xe and Ne arc lamps. For each 
night 100 bias and 100 well exposed dome flat-field images 
were taken (~30,000 ADU/pixel/image). Our program uses 
dome-flats, as opposed to sky-flats, as small scale pixel-to- 
pixel variations are a potentially important source of noise, 
requiring a large number of well exposed fiats to effectively 
remove it. Our observations were all specified to occur at 
the same pixel location, to help suppress flat-fielding errors. 
Through ESO GTC program 182.C-2018, we used 12 tran- 
sit light curves taken with narrowband imaging to estimate 
the guiding performance during a typical 4.5 hour long tran- 
sit observation. We found the guiding performance of GTC 
OSIRIS to be quite stable, with measured drifts of 1.4±1.0 
pixels (0.18±0.13") in the CCD X-direction and 1.0±0.6 pix- 
els (0.13±0.06") in the Y-direction over a 4 hour period. 



2.3 Reduction 

The bias frames and flat fields were combined using standard 
IRAF routines, and used to correct each image. Aperture ex- 
traction for both the target and reference spectra was done 
using a wide ~40 pixel aperture and background subtrac- 
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Figure 2. Extracted R500B and R500R grism spectra of XO-2A 
(blue & red) and the reference star XO-2B (black). Also overplot- 
ted (grey) is a R500B spectra of the sky background, increased 
in counts by a factor of 10 for illustration purposes. 



tion (see Fig. [2J . For both nights, the appropriate aperture 
size was chosen to minimise the out-of-transit light curve 
dispersion. The extracted spectra were then used to create 
differential spectrophotometric light curves, summing the 
spectra over wavelength regions of interest, and dividing the 
flux of the target star by that of the reference. With a low- 
resolution grism in conjunction with a wide slit and modest 
defocus, the spectra are inherently of low spectral resolution 
limiting bandwidths to ~30 A or larger (R=200), with 50 A 
bands the narrowest chosen for this study (R=120). 



3 ANALYSIS 

We modelled the transit light curves with the analytical 



transit models of Mandel & Agol (20021, fitting the inte- 
grated light curves for the central transit time, planet-to- 
star radius contrast, inclination, stellar density, and base- 
line flux using a Levenberg-Marquardt least-squares algo- 
rithm (Markwardt 2009). To account for the effects of limb- 
darkening on the transit light curve, we adopted the three 
parameter limb-darkening law, 



= 1 - C 2 (l -fjt)- C 3 (l - M 3/2 ) - C4(l 



M ) 



(1) 



calculating the coefficients following |Sing| ( |2010[ ). The er- 
rors on each datapoint were set to the standard deviation of 
the fit to the out-of-transit residuals. Following |Pont et al.| 
(2006), we assessed the levels of residual red noise by mod- 
elling the binned variance with a a 2 = cr„ /N + a 2 relation, 
where a m is the uncorrelated white noise component, N is 
the number points in the bin, and oy characterises the red 
noise. We subsequently rescaled the fitted parameter un- 
certainties to incorporate the measured values of ay in the 
white-light and spectral bin fits. 

In our fits, we determined the best-fit system parame- 
ters (inclination, central transit time, ect.) from the white- 
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Figure 3. The whitelight transit light curve using the R500B 
grism (black). Also overplotted is the fitted model light curve 
including parameterised systematic trends (red) and the residuals 
with error bars (black points bottom). 



light curve before measuring the differential radius varia- 
tions of the transmission spectra using wavelength bins. 



3.1 Broadband transit light curves and slit losses 

We summed up the R500B spectrum from 4066 to 7590 A to 
characterise the Feb. transit white-light curve (see Fig. [3). 
Summed over a wide spectral area, each of the white-light 
curve points has a very high total count rate, typically 3x 10 s 
counts per image, which gives a theoretical photon-limited 
precision of 7xl0 -5 . The observed transit light curve, how- 
ever, displays clear mmag-level systematic trends which vary 
on a ~2.5 hour timescale. We attribute these trends to dif- 
ferential slit losses between the target and reference star. 

The systematic slit-loss trends of the Feb. transit light 
curve exhibit a regular sinusoidal trend (see Fig. |4|. In dis- 
cussions with the GTC staff, the probable cause of the regu- 
lar and repetitive oscillation observed is the telescope guid- 
ing (A. Cabrera-Lavers priv. comm.). As the two objects 
are orientated following a specific slit-position angle and the 
guiding is performed with respect to the target, there are 
subtle slit losses due to the lack of rotator correction in the 
guiding. The guiding corrects by way of a linear displace- 
ment in the focal plane, and as no additional rotation is 
done, the effect in the non-target star will be translated as 
a small positional displacement that can be approximated 
as sine and cosine functions in the focal plane. With this 
interpretation, the regular trends in Fig. [4] are the result 
of differences in slit losses between the target and reference 
star, primarily due to guiding errors. 



3.2 Parameterised slit losses 

To model the systematic trends of the Feb. R500B light 
curve, we parameterised the slit losses with a combination 
of a sine and cosine function as a function of time, fitting for 
the sinusoidal amplitudes a, periods p, and phases 6. Upon 
searching for trends in other external parameters, such as 
CCD position and sky brightness, a linear full width at half 
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Figure 4. The Feb. R500B light curve with the transit removed 
illustrating the slit-losses (black) also overplotted with a param- 
eterised model (red). 



maximum (FWHM, S w ) trend and linear base trend in or- 
bital phase, (f>, were added as well. Shown in Fig. [4] the total 
time-FWHM-dependent baseline flux f(<f>, 5 W ) was described 
by a fit with nine free parameters, eight describing the sys- 
tematic trends and one measuring the differential system 
flux f , 



f(<f>,6 w ) = f X (ci0 + 1) x (oi sm(2-K(f)/p 1 + 0i) + l)x 

(oa cos(27T0/p 2 + e 2 ) + 1) x (c 2 S w + 1) (2) 

As with the sinusoidal trends, a FWHM dependance is also 
expected for slit losses, as wider Point Spread Functions 
(PSF) lead to wider wings and larger slit losses. 

The cosine term captured a higher harmonic of the over- 
all 2.5 hour modulation, as the fitted parameters were close 
to half the period, half the amplitude, and 90 degrees out 
of phase compared to the fitted sinusoidal terms. Further 
higher order (higher frequency) harmonics are evident in the 
residuals, which have an amplitude of ~2xl0 -4 with peri- 
ods all shorter than ~13 minutes. The sine and cosine terms 
effectively remove the hour-long timescale trends, relevant 
for measuring the planetary radius. 

With the parameterised slit losses, the standard devi- 
ation of the residuals is estimated to be 2.3xl0 -4 for the 
R500B white-light curve, which is dominated by the resid- 
ual slit-loss trends and a factor of three above the theoretical 
photon noise limit. We find a system inclination of 88.8±0.9 
deg and stellar density of 1.5±0.1 g cm -1 for the R500B 
transit (see Table [T]). These values match well with previous 
determinations, including the GTC photometric measure- 
ments of |Sing et al. ( |2011a l who found 87.62±0.51 deg and 
1.328±0.088 g cm" 1 . 

To measure the transmission spectra, we also sub- 
divided the R500B spectra into different sub-bins of ~50 A 
and wider. Parameterising the slit losses comes close to pro- 
viding a reliable transmission spectrum, which would be 
preferable as absolute transit depths would be measurable. 
However, we find that residual slit losses are still present at 
an unacceptably high level. The relevant radius precision one 
needs to observe atmospheric features is about one atmo- 



spheric scale height, estimated to be H—370 km (Sing et al 
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Table 1. Best fit parameters for the XO-2b R500B transit 



1 ';ii'Ai n n^r 1 ! 1 


Value 


Period, P< [days] 


2.61586178±0.00000075 


Mid-transit Time, [BJD] 


2455981.46036±0.00013 


radius, R p ;/R* 


0.1052±0.0011 


inclination, i [deg] 


88.8±0.9 


system scale, a/R* 


8.17±0.08 


impact parameter, 6=acosi/R* 


0.174±0.060 


stellar density, p* [g cm -3 ] 


1.507±0.1 


limb-darkening coeff., C2 


1.3036 


limb-darkening coeff., C3 


-0.8579 


limb-darkening coeff., C4 


0.3005 



Previously determined model dependent parameters 



R* [R ] 
R P i [Rjup] 
M* [M ] 
M p i [M Jup ] 



0.976±0.020t 
0.996±0.025* 
0.971±0.034* 
0.565±0.054* 



tfrom Sing ct al. (2011), fit along with the data from 
j2007| > and|Fernandez et al.| J2009| > 
iFromlFernandez et al.l J2009 i 



Burke ct al. 



2011a I. This target gives a required planet-to-star radius 
contrast precision of i//i?^=0. 00054 which is an absorption 
depth of A(7? pi /.R*) 2 =0.01%. With the adopted parameter- 
isation, the residual slit-loss trends and limiting precision 
is about twice this target value. Further observations with 
wider slits will reduce these observed slit-loss trends. With 
GTC, based on the observed pointing accuracy and PSF in 
good seeing conditions, 10" slits or wider should encompass 
the entire PSF and lead to minimal slit losses. 

The R500R observation (see Fig. ^ was best fit with a 
baseline flux and a linear trend with phase, a linear trend 
with X-position, plus a quadratic trend with FWHM. How- 
ever, the fits are of poor quality as the slit-loss trends are 
of much higher frequency, and are difficult to fit without 
including a large number of modes. 

3.3 Common-mode slit losses and differential 
transmission spectra 

As both the R500R and R500B lightcurves show multiple 
higher-order modes, requiring a prohibitory large number 
of free parameters, we also employed a common-mode tech- 
nique to remove the slit-loss trends. This technique relies on 
the similarities of the slit-losses across a spectrum, which can 
be characterised by the light curves themselves and removed 
individually for each spectral wavelength bin. For both 
the R500B and R500R, we find the overall systematics of 
the white-light curve dominate over wavelength-dependent 
trends by a factor of about six. Subtracting common-mode 
systematics are particularly useful when measuring trans- 
mission spectra, as the relative wavelength-dependent tran- 
sit depths are usually not effected. 

Empirically determining and removing the slit losses 
has an advantage over the parameterised method, as the 
higher order frequencies are naturally subtracted. However, 
absolute transit depths cannot be be obtained, and the slit 
losses do appear to have a slight wavelength dependence 
brought about through different flux illumination levels. 



1.005 



1.000 



=> 0.995 



S. 0.990 



0.985 



0.980 I ,_ 

-0.06 



-0.04 



-0.02 0.00 
Phase 



0.02 



0.04 



Figure 5. The whitelight transit light curve using the R500R 
grism. 

Higher count rates, at a particular wavelength, result in 
higher counts in the tails of the PSF, which spill over the 5" 
slit and lead to larger slit losses. 

We employed different techniques to remove the 
common-mode trends, including using the slit losses from 
the white-light curve (e.g. Fig. |4| as an auxiliary decor- 
relation parameter in the transit light curve fits (see Fig. 
[6j. We also used a differential-transit method, subtract- 
ing a reference-bin light curve from those of surrounding 
spectral bins, and fitting for the differential planetary radii 
ARpi/R s tar- In this method, the majority of the transit it- 
self is subtracted away, apart from the differences in limb- 
darkening and planetary radius, but the slit losses them- 
selves are also largely removed. Both methods produced sim- 
ilar results, though we quote the results and draw conclu- 
sions from the differential-transit fits, as we were able to fit 
the light curves with a minimum number of free parame- 
ters and avoided fitting for nuisance parameters describing 
the slit losses. The amplitude of the residual wavelength- 
dependent slit losses in the differential-transit fits were then 
estimated from ay for each wavelength bin, with the errors 
on the transmission spectra adjusted accordingly. 

The slit-loss dominated trends (e.g. Fig. [4| for the 
R500B and R500R grisms have standard deviations of 
62xl0~ 5 and 68xl0 -5 respectively, and thus have an am- 
plitude which is on the order of the size of the atmospheric 
spectral features we are searching for. The common-mode 
differential-transit technique, which removes a large portion 
of the slit losses, improves the limiting precision by about 
a factor of six, with typical values of the residual correlated 
noise measured to be oy ~ 10 x 10 -5 . The achieved preci- 
sion is about a factor of two better than the parameterised 
method, and sufficient to confidently measure and identify 
atmospheric features. 

We quantitatively assessed the spectral extent to which 
we could use the slit losses of a given reference bin by mea- 
suring the linear Pearson correlation coefficient between the 
out-of-transit portions of the reference and spectral-bin light 
curves. Well correlated slit losses produced correlation values 
of 0.8 and higher, giving red noise values near oy = 10 x 10~ 5 
in the differential-transit light curves. We found that spec- 
tral ranges were limited to no larger than ~500 A away 
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Figure 6. Three broadband transit light curves extracted from 
the R500B grism (black), correcting for slit losses using those 
from the white-light curve (as seen in Fig.[4j|, with an arbitrary 
offset for clarity. The effects of limb darkening are easily distin- 
guished from the bluest spectral bin (top, 4066 to 4745 A) com- 
pared to mid-spectral wavelengths (middle, 5849 to 6612 A) and 
the reddest bin (bottom, 7790 to 8582 A). The light curves have 
precisions of 400, 260, and 210 ppm with a cadence of 39 sec. (top 
to bottom respectively). 



from the centre of any particular reference bin, as correla- 
tion values drop significantly and the red noise increases. 
In the context of hot- Jupiter transmission spectra, this lim- 
ited our analysis to searches of expected specific atomic fea- 
tures including sodium, potassium, and H a . In our searches, 
potassium was unfortunately compromised by the large and 
variable telluric O2 lines (also see |Crossfield et al.|2012[ ), and 
the H a line did not show significant absorption. 

In the sodium search, we restricted our extracted spec- 
trum to the small region around the sodium feature. The 
wavelength ranges were 5450-6150 A for the R500B spec- 
trum and 5850-6150 A for the R500R spectrum. The effi- 
ciency of the CCD introduces significant errors in the R500R 
spectrum blue-ward of 5800 A, where the efficiency is less 
than half of that at the peak response. The light curves were 
fit adopting the stellar model limb darkening values as well 
as fixing the ephemeris, inclination, and stellar density to 
the best-fit values found previously. The fit included three 
free parameters, one parameter for the differential planet to 
star radius A R p i / R s tar and two parameters for the baseline 
flux, allowed to vary in time linearly. The two transits were 
fit individually, with both transits showing significant ab- 
sorption in the sodium band. The resulting combined spec- 
trum is shown in Figure [7J using bin sizes of 50 A. 

We tried excluding the reference star XO-2B from the 
differential-transit analysis, though found significantly worse 
results, indicating the reference star does help reduce the im- 
pact of slit losses in the photometric light curves. We also 
investigated the effect of changing the spectral binning on 
the resulting differential spectrum in increments of 12 A. In- 
cluding the variance of each point from this method had a 
negligible effect on the spectral shape or error bars, except 
in the sodium feature, where the absorption depth decreased 
if the band was moved more than 25 A away from the cen- 
tral wavelength of the feature, consistent with a positive 
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Figure 7. The differential spectrum around the sodium feature, 
using a weighted mean of both observations where the wave- 
lengths overlap (5850-6150 A). The X-axis error bars indicate 
the wavelength spectral bins, while the Y-axis error bars indicate 
the 1-tT uncertainty in differential radius. 



detection. In addition, wider bandwidths dilute the mea- 
sured sodium absorption depth, with the feature confined 
to bandwidths narrower than ~50 A. The red noise was 
seen vary by up to ~ 20 per cent depending on the partic- 
ular bins used, and we selected the binning which resulted 
in the lowest oy. We also found that changing the reference 
wavelengths by 50-100 A made no difference to the shape of 
the relative spectrum. The measured values for oy derived 
from the binning technique also correspond well with the 
observed dispersion of the transmission spectra, excluding 
the sodium line, giving further indications that the values 
are accurate. 

We also measured the differential absorption depth of 
just the sodium feature itself using the differential-transit 
method, which is similar to the procedures in Charbonncau 
et al. (20021. We measured the differential transit depth be- 
tween a 50 A band centred on the sodium doublet to a refer- 
ence composed of a blue (5812-5862 A) and a red (5917-5967 
A) band which bracket the sodium doublet. The differen- 
tial lightcurves were first normalised to the out-of-transit 
flux for the given wavelength ranges. This measurement is 
the most reliable way to compare the sodium absorption 
depths between the two transits taken with different grisms, 
as the same wavelength regions can be compared directly. 
We use this measurement to assess the significance of the 
sodium feature. We found significant and consistent absorp- 
tion depths between the two transits with different grating 
settings, measuring absorption depths of 47 ± 11 x 10~ 5 for 
the R500B observation and 49 ± 17 x 10~ 5 for the R500R 
observation (see Figure [8|. The combined relative absorp- 
tion depth of the Na I feature using the mean spectrum is 
47 ± 9 x 10 -5 , which is consistent with the spectral results. 



4 DISCUSSION 

4.1 Atmospheric model fits 

The spectrum appears flat around the sodium feature, sug- 
gesting that the doublet must be unresolved in the data, and 
narrower than 50 A. Such a narrow feature could be indica- 
tive of haze or clouds, which obscure the lower regions of the 
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Figure 8. Differential transit light curves of a 50 A band cen- 
tred on the sodium doublet to that of a reference composed of a 
blue (5812-5862 A) and a red (5917-5967 A) band which bracket 
the sodium doublet. Both transits using the R500B (black circles) 
and R500R (red squares) are plotted, binned by 4 and 14 points 
respectively along with the best-fit models (black and red re- 
spectively). The R500R light curve has stronger differential limb- 
darkening (producing a deeper, more rounded transit) due to the 
weaker spectral responce in the blue reference band. 



atmosphere by scattering signatures, as seen in HD 189733b 
(|Pont et al.||2008| |Lecavelier des Etangs et aT1|2008| |Desert| 



et al.|2009||2011||Sing et al.|2011b||Huitson et al.|2012[ ). Al- 

ternatively, in clear atmospheres dominated by alkali metal 
absorption, flat spectral regions can be observed where the 
abundance drops suddenly as altitude increases which coun- 
teracts the increase in cross-section toward the line core. 
This behaviour is seen in HD 209458b ( |Sing et al.1|2008} 
Vidal-Madjar et al.||20lfb l. Abundance drops could be at- 
tributed to processes such as ionisation, or condensation if 
the temperature profile crosses the condensation curve. We 
fit the spectrum of XO-2b with different models to assess 
the significance of the observed lack of line wings. 

We first used a model transmission spectrum developed 
specifically for XO-2b from |Sing et a l. (201 la]) which is calcu- 
lated using the hot- Jupiter model atmospheric calculations 
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Figure 9. The differential transmission spectrum from the 
R500R and R500B spectra around the sodium feature (black cir- 
cles). Also plotted (blue solid line) is a model fit from [Fortney| 
|et al.| l |2010| ) , assuming a haze / cloud free atmosphere dominated 
in the optical by alkali metal absorption. 



freely produced a temperature of 390 ± 120 K (x 2 = 31.6 for 
12 DOF) which is clearly too cold, showing that the model 
has to flatten the line wings to fit the data. Neither model 
provides a satisfactory fit. 

For an alternative model, we assumed that the sodium 
feature is unresolved at a resolution of 50 A with the line 
wings either obscured by scattering from clouds or hazes, 
or flat due to a sudden abundance change. We choose a 
width of 12 A for the unresolved feature, which must be less 
than 50 A, and then artificially set the regions of the high- 
resolution model outside this wavelength range to be flat 
with a constant radius. The flat regions of the model were 
set to the zero reference level and fit with a free parameter. 
In this model, all regions outside the 12 A feature have the 
same radius, which is allowed to move up and down relative 
to the sodium feature. This is equivalent to adjusting the 
height of the cloud deck/haze layer, or adjusting the sodium 
abundance at the lower altitudes probed by the sodium line 
wings. Additionally, the Zo parameter is still free, allowing 
the model spectral profile to freely move up and down as a 
whole. The resulting fits are the same when this procedure 
is employed on both the analytic isothermal models from 



of Fortney et al. (20101 and Shabram et al. (20111. These Huitson et al. ( 2012 \ and the numerical Fortney et al. ( 2010 1 



models include a self-consistent treatment of radiative trans- 
fer and chemical equilibrium of neutral species, however the 
formation of hazes and photo-ionisation are both neglected. 
Broad sodium and potassium line wings are seen, each over 
1000 A wide. The model was binned to the transmission 
spectral resolution of 50 A and fitted to the data with the 
reference level zo as a free parameter, as shown in Figure [9] 
The best fit gives X 2 = 37.6 for 13 DOF. 

We then modelled the feature using an analytic isother- 
mal model transmission spectrum for the sodium feature as 
described in|Huitson et aL](|2012[), b ased on the formalism of 



models, as the narrow sodium feature is only detected in one 
50 A bin, provid ing no temperature information. The fit is 
The best fit gives a = 19.2 for 12 



10 



Lecavelier des Etangs et al. (2008) and Vidal-Madjar et al 



( 2011b I, with temperature and the reference level, zq as free 



parameters. Fixing the isothermal model temperature to an 
average temperature of 800 K produced a fit very similar 
to the full model, with X 2 = 38.3 for 13 DOF. Adding pres- 
sure broadening to this model produced no noticeable differ- 
ence in the fits. Allowing the temperature parameter to vary 



shown in Figure 
DOF, which is a better fit than the standard model at the 
3.5 — 4a level. 

Changing the width of the unresolved feature, wider or 
narrower than 12 A did not change the best fit, since the 
free parameter for the altitude of the line wings adjusted to 
compensate. Additionally, we tried adding a Mie or Rayleigh 
scattering cross-section to the sodium line cross-section in 
the isothermal models, to try to distinguish between the 
two cases of an obscuring haze or a sudden abundance de- 
crease from photo-ionisation or condensation. We find that 
we cannot differentiate with any significance between these 
cases, with the resulting \ 2 values only la different. 

Both condensation and ionisation have been inferred 
in hot Jupiter atmospheres ( |Sing et al.||2008| |Vidal-Madjar| 



et al. 2011b I. Although the global average theoretical ter- 
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Figure 10. Fit assuming an unresolved feature narrower than 
50 A wide, either due to a cloud or haze layer obscuring the line 
wings, an abundance drop with increasing altitude. Fitted data 
points arc shown in black with circles and joined with dotted 
lines, and the model is shown with a blue solid line. 



minator T-P profile does not cross the condensation curve 
of Na into Na2S, it is possible for temperatures to be low 
enough for this condensation on the planet's night side, thus 
causing a decrease in Na I abundance in the upper atmo- 
sphere. 



5 CONCLUSION 

We detect absorption from the sodium doublet in the at- 
mosphere of XO-2b at the 5.2a level using multi-object 
long-slit spectroscopy. This makes XO-2b the first extra- 
solar planet where both potassium and sodium have been 
detected in its atmosphere, features which have both been 
long predicted to be present in hot Jupiters. The detec- 
tion helps demonstrate the feasibility of using ground-based 
multi-object spectroscopy for detecting atmospheric features 
in transiting extrasolar planets. While our analysis was lim- 
ited to measuring differential transmission spectra, due to 
slit losses, the precisions achieved on absolute transit depths 
(200 ppm) would already be sufficient to detect features on 
nearly a dozen exoplanets with larger scale heights, such 
as Wasp-17b, and adopting wider slits should substantially 
increase the photometric precision. 

Of the well studied planets so far, an emerging pic- 
ture is far from simple, with observations of clouds or hazes 
and large abundance changes throughout the atmospheres 
detected in the wide class of hot- Jupiter exoplanets. Fur- 
ther observations of XO-2b will be able to confirm and help 
distinguish the cause of the detected lack of sodium line 
wings, as a full broad-band transmission spectra could pick 
up Rayleigh scattering from haze, for instance. Addition- 
ally, these low-resolution observations would highly bene- 
fit from complementary high-resolution transmission spec- 
troscopy around the sodium and potassium features, as the 
measured alkali metal absorption line profiles can place con- 
straints on temperatures and abundances. 
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